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Intermediate Scattering Functions (ISF's) are measured for colloidal hard sphere sys- 
tems using both Dynamic Light Scattering (DLS) and X-ray Photon Correlation 
Spectroscopy (XPCS). We compare the techniques, and discuss the advantages and 
disadvantages of each. Both techniques agree in the overlapping range of scatter- 
ing vectors. We investigate the scaling behaviour found by Segre and Pusey^ but 
challenged by Lurio et alP. We observe a scaling behaviour over several decades in 
time but not in the long time regime. Moreover, we do not observe long time diffu- 
sive regimes at scattering vectors away from the peak of the structure factor and so 
question the existence of a long time diffusion coefficients at these scattering vectors. 
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I. INTRODUCTION 



Over the past decade, X-ray Photon Correlation Spectroscopy (XPCS) has emerged as a 
valuable tool for the study of dynamics in soft matter^!. The technique has two principle 
advantages over the well established technique of Dynamic Light Scattering (DLS): the 
first is that the shorter wavelengths allow for measurements at larger scattering vectors 
(smaller size scales); the second is that X-ray measurements do not suffer from multiple 
scattering. These advantages are balanced against two major disadvantages: X-rays have 
lower coherence than (laser) light sources (leading to reduced correlation intercepts); and 
that samples can be easily damaged by intense X-rays. 

A number of authors have applied XPCS to study the dynamics of suspensions of col- 
loidal particles^EJ] This builds on extensive studies using light scattering technique d 12 * 13 !, 
and, more recently, microscopy techniques"^. The most widely studied system is that of 
particles which behave as hard spheres". 

The fundamental quantity measured in scattering experiments is the Intermediate Scat- 
tering Function (ISF) which describes the dynamics of the particle number density. For 
suspensions at moderate concentrations, the dynamics are often divided into three regimes: 
a short time diffusive regime, where particles diffuse within their neighbour cages; a non- 
diffusive crossover regime or plateau, where the interactions between a particle and its 
neighbour cage are most clearly exposed; and a long time diffusion regime where the par- 
ticles have escaped their neighbour cages. This long time regime is particularly difficult to 
access experimentally, and ambiguities in interpretation can arise depending on the method 
of analysis used. 

To date there have been no direct comparisons between the two techniques on concen- 
trated suspensions of colloidal particles. There are however instances where XPCS and 
DLS experiments have not been in agreement. Specifically, DLS experiments carried out on 
sterically stabilized hard spheres^, found that the short and long time diffusion coefficients 
(D s (q) and D L (q) respectively) both scale with the structure factor - in other words the ratio 
D s (q)/DL(q) is independent of q. However, this finding has since been questioned by Lurio 
et al.™ who used XPCS to study a system of charge-stabilized polystyrene latex spheres in 
glycerol. In this work they showed that the measured structure factors for this system were 
consistent with those of hard spheres - however, when they investigated the dynamics of 
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this system, they did not observe the long time scaling. There are several possible expla- 
nations for this discrepancy: i) there is a fundamental difference between XPCS and DLS 
which results in different measured behaviours; ii) X-ray damage contributes to the mea- 
sured results in XPCS; iii) the charge- stabilized pseudo-hard-sphere system is not equivalent 
to the steric hard sphere system; and/or iv) the differences may arise from ambiguities in 
the determination of the long time diffusion coefficient! 18 ' 19 ^. 

In this paper we compare the two techniques by examining sterically stabilized hard 
spheres using both DLS and XPCS. By using the same samples for both techniques we 
can eliminate sample variation and effects of sample preparation. We conduct analyses 
similar to those in Lurio et alP and Segre and Pusey^, and compare the results of the two 
methods. A previous study^ investigated the comparison between DLS and XPCS but only 
at a low volume fraction, 0.164. In that work the overlapping range of scattering vectors 
did not include the peak of the structure factor or the minimum of the form factor, where 
contributions from multiple scattering are maximal. 



II. THEORY 

DLS and XPCS are now standard methods for measuring the dynamics of colloidal sus- 
pensions under certain conditions and have been reviewed in numerous article d 4 ' 7 ' 12 ' 19 ™ . In 
this section, we summarize some important aspects. The basic quantity measured by both 
techniques is the (normalized) time-averaged autocorrelation function (ICF) of the intensity, 
I(q,t), of the scattered light: 

where q is the magnitude of the scattering vector, r is the delay time. For ergodic systems 
the ICF factors according to the Siegert relationship: 

g {2 \q,r) = l + c\f(q,r)\ 2 (2) 

where 



and 



F (Q, r) = i £ (exp [iqir-ir) - r k (0))}) (4) 

j,k=l 



is the intermediate scattering function - the auto-correlation function of the q th spatial 
Fourier component of the particle number density fluctuations. The polydispersity of the 
colloidal suspension, i.e. the (small) relative spread in particle radii (see below), is not 
explicitly taken into account. In Eq. ([2]), c is an experimental constant determined by the 
ratio of the coherence area to the detector area. Note also that F(q, 0) = S(q) is the static 
structure factor. 

For diffusive density fluctuations we have, 

f(q,r) = eM-q 2 D(q)r) (5) 

More generally one calculates the time-dependent quant ity^^ 

= (6) 

The latter describes the evolution from short time processes, expressing the particles' (dif- 
fusive) response to those hydrodynamic modes in the suspending liquid that propagate 
instantaneously on the experimental time scale r > 10 _6 s, to the long-time processes asso- 
ciated with structural rearrangement or diffusion over larger distances. These processes are 
characterised by the short and long time collective diffusion coefficients respectively: 

D s {q) = lim [D{q,r)} and D L {q) = lim [D{q, r)] (7) 

Where 77 is the lower limit of the experimental time window. It is now well established 
that D s (q) shows a minimum at the position, q m , of the maximum in the static structure 
factor, S(qj^. This is the familiar de-Gennes narrowing - the initial diffusive decay of density 
fluctuations is slowest at gj 22 ' 23 1 . 



III. METHODS 

A. Sample Preparation 

The particles used in this work consist of a co-polymer core of methymethacrylate (MMA) 
and Trifluorethylmetacrylate (TFEMA). The co-polymer TFEMA is added to enable refrac- 
tive index matching in a single solvent (cis-decalin), and the suspension turbidity can be 
controlled by changing the temperature. For these samples the index matching is suffi- 
ciently good that multiple scattering is nearly zero in the range of scattering angles studied 
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here. The particles are sterically stabilised by a thin coating of poly-12-hydroxystearic acid 
(PHSA), approximately 10 nm thick. The particles used here are designated XL63, and 
have a radius R = 185nm and a polydispersity of 8 — 9%pH2D All times are presented in 
units of the Brownian time, tb = R 2 /QD Q = 0.013s where D Q is the diffusion constant for 
freely diffusing particles. The hard sphere-like interaction between these particles in these 
suspensions is confirmed by the identification of the equilibrium phase boundaries, i.e. the 
freezing and melting points. 

Samples for DLS were prepared in 8mm path length, sealable scattering glass cells. Sam- 
ples for XPCS were transferred to quartz capillaries of 1.5mm diameter (Wolfgang Muller 
Glas Technik, Berlin). A ball bearing was placed into the cell to facilitate the tumbling of 
the particle suspension to ensure good mixing. The cells were sealed with araldite, ensuring 
evaporation was insignificant over timescales of months. For some test measurements (see 
below) DLS measurements were also carried out in the X-ray capillaries prior to or following 
the X-ray measurements. 

An ALV-6010 spectrometer is used for the DLS measurements. A HeNe laser of wave- 
length 632.8 nm illuminates the sample, and an Avalanche Photon Detector (APD), located 
on a goniometer, measures the scattered intensity at a specific scattering angle 6. The acces- 
sible angular range is 15 — 150°, corresponding to a q range of ~ 3.8 — 28 / um~ 1 , encompassing 
the peak in the structure factor at q ~ 19/xm -1 . 

The XPCS experiments were performed at the IDlOa station of the Troika beamline at the 
European Synchrotron Radiation Facility (ESRF) in Grenoble, France. A third generation 
synchrotron, combined with an undulator insertion device, provides a partially coherent X- 
ray beam which has sufficient brilliance for XPCS experiments. The standard "coherent" 
set-up is based on two factors: (1) a strong collimation of the beam to match the transverse 
coherence length £t of the TROIKA undulator source, £t = XL/s, where s is the source size 
and L the sample-to-source distance; and (2) wavelength selection using a standard Si(lll) 
monochromator that provides the longitudinal coherence length £^ = A 2 /AA. At the ID10A 
beamline this routinely yields a 10x10 /im 2 beamspot at the sample with a 1/xm longitudinal 
coherence length at 8 keV (A = 1.555 A), well inside the limit for path length differences 
in the present range of scattering vectors cp. Focusing optics (compound refractive lenses) 
can be inserted into the beam to optimise the intensity whenever the coherence conditions 
can be relaxed. In this configuration the "coherent" intensity on the sample can exceed 
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10 10 photons/s. To match the speckle size a OD detector (scintillation counter) is placed 
behind a 100 \xm aperture at a distance D = 2240mm behind the sample. Measurements of 
the intensity fluctuations at small angles 9 = d/D are achieved by translating the detector 
by a distance d perpendicular to the incident beam. A digital autocorrelator, connected 
to the detector, allows the direct measurement of the ICF at a specific scattering vector of 
amplitude q in the range ~ 10 — 100/im" 1 . Capillaries were mounted in an evacuated sample 
chamber with temperature control (T = 21°C). 

IV. RESULTS 

Fig. [T] shows a comparison between the two techniques for a dilute suspension of hard 
sphere particles for three values of qR. Fig. [IJi shows the ICF (Eq. [2]), which demonstrates 
the lower intercept of XPCS (symbols) relative to DLS (lines) - this is a consequence of 
the lower coherence of the X-rays relative to visible light. Fig. shows |/(g,r)| 2 , after 
normalization by the intercept, demonstrating that the two techniques give very similar 
results, despite the difference in intercept. 

Fig. [2] shows the ISFs for both techniques. As can be seen, the reduced intercept of the 
XPCS measurements now becomes important, and leads to the non-zero baseline for two of 
the scans. This is exacerbated by poor statistics, due to the fact that at low volume fractions 
the scattering is very low. This could be improved by accumulating for longer, however in 
the present case this was not attempted to limit the possibility of beam damage. Despite 
this difference, the results of the initial decay are still consistent. 

DLS measurements were typically carried out in large sample cells with ~ 8mm path 
lengths, whereas X-ray measurements must be carried out in small path length (~ 1.5mm) 
X-ray capillaries. To determine whether the use of capillaries causes any change in the 
dynamics (eg due to shear alignment during loading), the inset in fig. [2] shows a comparison 
between DLS measurements at two volume fractions in standard DLS cells and in X-ray 
capillaries. The agreement seen here demonstrates that any difference observed between 
DLS and XPCS results is not due to different preparation methods or sample cells. 

Turning now to higher volume fractions, the measured ISFs from both techniques are 
shown in fig. [3] for a range of volume fractions and scattering vectors. Fig. [3] shows that, in 
general, there is reasonable agreement between the two techniques over the range of volume 
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Figure 1. XPCS (open symbols) and DLS (lines) measurements of the intensity autocorrelation 
functions for a dilute suspension (</> ~ 0.01) of XL63 at the qR values indicated; (a) un-normalized: 
(g( 2 \q,T) — 1) and; (b) normalized: (g( 2 \q,r) — l)/c. 




Figure 2. Calculated ISFs from XPCS (open symbols) and DLS (lines) measurements on dilute 
suspensions at the qR values indicated. The inset shows DLS measurements carried out in DLS 
cells (lines) and X-ray capillaries (open symbols) at the volume fractions indicated at a qR value 
of 3.48 (near the peak of the structure factor). 
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Figure 3. ISF from XPCS (open symbols) and DLS (lines) at the volume fractions and qR values 
indicated. 

fractions shown, though deviations can be seen, particularly where the statistics are such 
that XPCS does not reach a zero baseline. 

In order to compare the two techniques in more detail, we fit stretched exponential 
functions with the form f(q, r) = exp[— (t/t£)^)] by minimizing the Chi-square value using 
a form of non-linear least squares (Levenberg-Marquardt algorithm) with Igor Pro software 
from wavemetrics. DLS data are fit out to delay times t/tb = 20. Due to the increased 
noise, XPCS data were fit over a smaller range, out to t/tb = 5. 

The results of this analysis are shown in fig. [4] for both DLS and XPCS data. Clearly 
there is agreement between the two techniques at volume fractions of 0.318 and 0.475. At 
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Figure 4. (a) Characteristic time and (b) exponent /3 of stretched exponential fits to the ISF, 
shown as functions of qR, at the volume fractions indicated. XPCS (filled symbols) and DLS (open 
symbols connected by lines). Error bars are estimated from the fits. For the characteristic times, 
the error in the Brownian time is also taken into account. Errors bars for DLS are smaller than 
symbols. Lines are drawn to guide the eye. 



a volume fraction of 0.421, XPCS gives characteristic times slightly larger than DLS and 
values slightly lower than DLS. However, the trends are clearly the same. The XPCS 
data, although it is subject to higher statistical noise, gives reasonable results out to higher 
qR values. Most importantly, the two techniques show the same variation with q. Having 
established that XPCS and DLS data give consistent results we now turn to the question of 
scaling. 

The above results show that the differences in scaling observed by Lurio et alP and Segre 
and Pusey^ cannot be due to differences between the two techniques. Further, although the 
samples can become damaged if exposed to the X-ray beam for too long, these results show 
that, with care, damage can be avoided, suggesting that this is not an explanation for the 
differences observed in previous measurements. However, the possibility of beam damage 
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Figure 5. Long time (top) and short time (bottom) behaviour of the ISF determined using DLS 
at (j> = 0.456 for several qR values as indicated. The representation is the one used in Segre et 
affl The linear behaviour at short time allows the identification of a short time diffusion coefficient 
D s (q). 

provides an upper limit on how many runs can be made on each sample. By contrast, for 
DLS, an arbitrarily large number of runs can be made to achieve the desired statistical 
reliability. For the data studied here 50 runs of 1000s were routinely made. For this reason, 
having established the equivalence of the two techniques, the analysis below will be limited 
to the DLS data. 

Fig. [5] plots ln(ISF) vs q ,2 r/r^ at a volume fraction of 0.456 for a range of scattering 
vectors for long (top) and short times (bottom). Note that the range of delay times shown 
is similar to those presented respectively in fig. [iji and fig. [3^ of ref (1), and the noise on 
the DLS data is very low, demonstrating that our data are comparable with those of Segre 
et al. 

Segre et alP and then Lurio et alP estimated short and long time diffusion coefficients 
by fitting straight lines to the initial and final decays of ln(ISF) vs q 2 r, as shown in Fig. 
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[5| Such analyses inherently assume, from the outset, that the fastest and slowest detected 
density fluctuations are diffusive. The validity of this assumption will be discussed below. 
The full identification of the short and long time diffusive regimes requires the observation of 
plateaux in the time-dependent diffusion coefficient D(q,r) (eq. 5) at short and long times 
respectively. To avoid duplicating unnecessary data, we present directly D(q,r)/D s (q), 
where D s (q) were obtained by identifying a plateau in D(q,r) at short times. The results 
are shown in fig. [6] for several scattering vectors at 6 volume fractions. The plateaux at 
short times suggest that the assumption of the existence of a short time diffusion coefficient 
is reasonable. However, this is not the case at long times, where there is no clear plateau in 
most cases, and thus the long time diffusion coefficient cannot be unambiguously determined. 

Despite the noise, our results are in qualitative agreement with fig 2a of ref (1). A scaling 
of D(q, r)/D s (q) is reasonably observed at volume fraction <fi > 0.456 for qR values spanning 
from 1.85 up to 4.91, even though noise is present at <ft — 0.456. As the volume fraction is 
increased the q dependence appears to decrease and the curves begin to converge. This trend 
is clearest for the volume fractions of 0.456, 0.498 and 0.531. At 0.456 the curves converge, 
within the noise, except for the lowest two qR values of 0.94 and 1.40. At a volume fraction 
of 0.498, the data for qR = 1.85 begins to deviate too, and at 0.531, in the metastable regime 
(ie where crystallization will occur at sufficiently long times), the deviation of the lowest 
three data sets becomes clearer. Interestingly however, while the deviation of the data for 
the lowest scattering vectors becomes more pronounced as volume fraction increases, the 
scaling for the scattering vectors around the peak becomes better. This scaling fails at long 
times, as highlighted in the log-log representation shown in the insets of fig. [6]for the highest 
volume fractions. By contrast with the results presented here, Segre et al. clearly observed a 
plateau at long times at <fi = 0.456 for all qR values spanning from 1.0 to 3.9 while we observe 
such plateaux only at qR values near the peak of the structure factor. Thus we question the 
existence of a long-time diffusive regime, and so a long-time diffusion coefficient D^q), in 
our system. 

The numerical derivative required to calculate D(q,r), introduces noise and makes it 
difficult to identify a long time diffusive regime, where such a regime exists. An alternative 
analysis which is not subject to this noise, and more clearly exposes the diffusive regimes, 
is to calculate the width functional, defined by analogy with the mean square displacement 
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Figure 6. Time-dependent diffusion coefficient D(q,r)/ D s (q), from DLS experiments, normaiised 
by the short-time diffusion coefficient at qR values and volume fractions indicated. Insets show 
same data on a double-log scale. 

as: 

w(q,r) = -\n[f(q,r)]/q 2 (8) 

which is shown as a function of delay time in fig. [7| The insets show the width function at 
q m . In this representation diffusive regimes are identified where w(q,t) grows linearly with 
delay time. This determination of diffusion does not require numerical differentiation of the 
data (in contrast to fig. [6]), and is insensitive to the time scale on which the ISF is plotted (in 
contrast to fig. [5J. Clearly, this representation shows unequivocally that short time diffusive 
regimes are observed over several decades for all volume fraction and qR values probed. 
However fig. [7] confirms that difficulties occur in identifying linear long-time behaviour at 
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some q values. In other words, the ISF and so the width function attain the experimental 
noise floor before the anticipated long-time diffusive regime is observed. For the data shown 
here, long-time diffusive behaviour can only be discerned at q m for > 0.498, and this is 
shown in the insets - the long time diffusive regime is not attained for other q values or for 
lower volume fractions. This is consistent with recent experiments which showed that the 
mean-squared distance particles must traverse, in order for density fluctuations to forget the 
effects of packing constraints, is smallest at qdr^. A consequence of this is that we are unable 
to directly conclude whether or not the ratio D S /D L is independent of q, as found by Segre 
et alP, a result which was disputed by Lurio et alP. Returning to the scaling behaviour, as 
(ft is increased, D(q, r)/D s (q) seems to scale better in the first decades even though the last 
decades shows an increasing q-dependence once > 0.498. To quantify this observation, we 
choose three fixed delay times, and plot in fig. [8] the standard deviation of the spread of 
D(q,r)/ D s (q) values at these times as a function of volume fraction. The scaling is clearly 
observed via the convergence of this quantity as volume fraction is increased. Of course the 
decay of the ISF shifts to longer delay times as volume fraction increases, so this analysis is 
limited. However it clearly demonstrates the scaling as volume fraction increases. Keeping 
in mind that the determination of the plateau at long times is ambiguous, particularly at 
low volume fractions, fig. [9] shows the qR dependence of the inverse of the short time 
collective diffusion coefficient. Clearly D /D s scales with the structure factor in agreement 
with previous results eg Segre and Pusejffl We stated above that the determination of the 
long time diffusion coefficient is problematic. In order to remove the ambiguity associated 
with the measurement of a long time plateau, we simply determine the apparent diffusion 
coefficients D(q,r x ) at several delay times t x . The ratio D(q,r x )/D s (q) is then plotted in 



fig. [10] for several volume fractions. At the lowest volume fraction measured there is some 
hint of scaling only at the early times but clearly not at long times. As the volume fraction 
increases any hint of scaling disappears, and the ratio never demonstrates q-independence 



at any delay time. Finally, Fig. 11 shows D(q,r)/D s (q) near the peak in the structure 
factor for several volume fractions. This graph confirms that at the peak it is just about 
possible to define a diffusive regime at long times (see also fig. [7]f). Of more interest is the 
volume fraction dependence of the rate at which the behaviour deviates from diffusive. At 
t/t b = 10° (after three decades in time) D(q,r)/D s (q) has dropped by almost an order of 
magnitude at a volume fraction of 0.549, whereas D(q, r)/D s (q) remains approximately 1 at 
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Figure 7. Width function as function of delay time at qR values and volume fractions indicated. 
Lines and dashed lines are linear fits of log[w(q, r)] = log(D s {q)) + log(r) at, respectively, short 
and long times. Inset shows the width function at the qR values nearest to the peak of structure 
factor for the corresponding volume fraction (scale as main panel). 

a volume fraction of 0.213. So, as the volume fraction increases, the system deviates more 
quickly from diffusive behaviour, even though the dynamics are slower. 



V. DISCUSSION 

The first part of the results compared DLS and XPCS results on identical samples. As 
far as we are aware, this is the first comprehensive comparison between DLS and XPCS on 
concentrated hard sphere samples. The results demonstrate a number of things. Within 
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Figure 8. Standard deviation of D(q, t/tb)/ D s (q), from DLS experiments, at fixed values of t/tb 
and for qR values in the range [2.29,4.91]. 
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Figure 9. Inverse of the short time diffusion coefficient expressed in term of the free diffusion 
coefficient Dq at volume fractions indicated. Lines are drawn to guide the eye. 

statistical errors, the ISFs determined from DLS and XPCS agree very well. The main 
limitation of XPCS is that good sample statistics are much harder to obtain, partly due to 
the lower coherence, and partly due to beam damage, which limits how long statistics can 
be accumulated for on a particular sample. This is particularly true for the higher volume 
fractions, where restricted diffusion means individual particles are exposed to the X-ray 
beam for long periods of time. In order to improve the statistics at high volume fractions, 
it would be desirable to have the samples mounted on a translation stage so that a series 
of short measurements can be made at different positions in the sample to limit sample 
damage. Such studies are currently being planned. 

Despite this, the results confirm that careful use of XPCS can provide data of similar 
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Figure 10. Ratio of D(q, T x )/D s (q) as a function of qR at volume fraction <\> equal to 0.475 (circles), 
0.498 (squares), 0.519 (triangles), 0.531 (diamonds), 0.540 (stars) and 0.549 (crosses). Each graph 
corresponds to a specific r x value as indicated. The bottom panel shows D(q,r x )/D s (q) at the 
longest time before noise starts to dominate. Note that this time increases with volume fraction. 
Within the noise there is no real indication of scaling between short and long time behaviour. 

quality to DLS. There are a number of situations where XPCS would be preferable to DLS: 
first, where particles are smaller than R ~ 120nm, where the structure factor peak cannot 
be accessed using DLS; second, for samples where multiple scattering is significant; and 
finally, where dynamics need to be studied at higher q values as shown in fig. [4] (though not 
explored here, even higher values are possible using XPCS). 

We turn now to the results of the scaling analysis. First, we find that D s (q) scales with 
S(q), an effect knows as DeGennes narrowing. Second, we find that the long time diffusion 
coefficient cannot be defined away from the peak for our system, so the scaling of the ratio 
D S /D L observed by Segre and Pusejff, but not by Lurio et alP, cannot be confirmed. Third, 
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Figure 11. Time-dependent diffusion coefficient, normalised to the short-time diffusion coefficient, 
as a function of delay time for several volume fractions at qR = 3.57 (around the peak of the 
structure factor for higher volume fractions). 



we do observe a scaling behaviour of the time-dependent diffusion coefficient at high volume 
fraction. This result agrees with Segre and Pusey but is not observed by Lurio et al. 

As the system under study is a hard sphere system very similar to the system used by 
Segre and Pusey, (though with higher polydispersity) the difference in the results is puzzling. 
One possible explanation is in the way the long time diffusion coefficient is determined. Our 
results have shown that the determination of this quantity is problematic, and according to 
our analysis, it cannot be defined away from the structure factor peak. Inspection of Fig. 



11, for example, shows that a clear plateau can be identified at early times (D s ), but that 
a plateau is much less well defined at long times. It is possible that the differences between 
our results and those of Segre and Pusey are related to the determination of this property. 

The identification of a long-time diffusive regime for the collective dynamics (through 
the width function, fig. [7]) is clearly observed around q m but not at other q values. So the 
following question can be asked: does a long-time diffusive regime exist for all ql We consider 
the quantity (Ar 2 (r^(g))), which corresponds to the mean squared distance a particle has to 
move for the number density fluctuations to forget excluded volume effects - where (Ar 2 ) and 
T m(l) are respectively the time-dependent mean square displacement and the delay time at 
the crossover between the fast and slow processes. This quantity was previously introduced 
and measured for the same particle suspensions^. At 0/ = 0.498 and q m , the quantity 
(Ar 2 (r^(g m ))) is smaller by a factor of ~ 5 compared to (Ar 2 (r^(g))) at qR=1.0. This 
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quantifies the difficulty in identifying the long-time diffusive behaviour from the coherent 
ISF at q vectors away from the peak of the structure factor, simply because particles are 
unable to move the distances necessary for the number density fluctuations to forget packing 
constraint effects. 

Interestingly, this question has been independently addressed in a recent paper by 
Holmqvist and Nagele^ for charge stabilized colloids. These authors do observe the long 
time scaling behaviour found by Segre and Pusey^ for scattering vectors around the peak 
of the structure factor. One difference between their system and the one studied here is 
that due to the smaller size, lower viscosity and perhaps the nature of the interactions, the 
timescales are very much faster - the ISF is complete within « 30 ms, whereas for the data 
presented here the decays at the highest volume fractions are of the order of 10 s, almost 
1000 times longer. In our experiments 50 x 1000 s measurements were made, but we were 
still unable to observe long time diffusive behaviour except at the peak of the structure 
factor. So it may be that for the hard sphere system used here, the determination of a true 
long time limit away from the structure factor peak is not possible with any reasonable 
experimental time. Nevertheless, an analysis similar to our figure 7 for the data shown in 
Holmqvist and Nagele^ would provide a more unambiguous measure of the determination 
of the long time diffusion coefficient, and would shed more light on the observed differences. 



More generally, the results summarized in fig. 11 show that the evolution of the time- 
dependent diffusion coefficient changes dramatically with volume fraction. At low volume 
fractions the transition between the two takes a long time. However as volume fraction 
increases the short and long time processes quickly separate, as more and more particles 
become trapped in neighbour cages. The magnitude of this separation increases, and occurs 
at shorter times, as the dynamics become slower. 



VI. CONCLUSIONS 

We present an extensive study of concentrated colloidal hard sphere suspensions using 
both DLS and XPCS, as functions of volume fraction and scattering vector. A scaling be- 
haviour, found by Segre et alP, is observed over several decades in time but not in the long 
time regime where significant q-dependence is observed for higher volume fractions. While 
the short-time diffusion regimes clearly exist at all volume fractions and scattering vectors, 
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long-time diffusive regimes are only observed near the peak of the structure factor at high 
volume fractions. Thus the existence of a collective long-time diffusive regime is rather ques- 
tionable, even though an apparent linear regime is observed in the common representation 
ln(ISF) vs q 2 r used by both Lurio et al 2 . and Segre and PuseyM By interpreting the quan- 
tity (Ar 2 (r^ l (g))), measured by van Megen et al.^, we explain the absence of a long-time 
diffusive regime by simply arguing that particles are unable to move the distances necessary 
for the number density fluctuations to forget packing constraint effects during accessible 
time scales. 
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